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SUMMARY 

The lack of a general theory for predicting turbulent flows has resulted in the 
development of various empirical techniques applicable to specific classes of these 
flows. One class of flows of considerable interest for many years, because of the 
various engineering applications, is designated as free turbulent shear flows. A gen- 
eralized eddy-viscosity approach has been successfully applied to these flows and is 
reported. Results presented herein for the test cases selected for evaluation by the 
Data Selection Committee of the NASA Working Conference on Free Turbulent Shear 
Flows show that predictions were obtained which are adequate for most engineering 
applications. 

Because of the importance of starting computations from the injection station 
where experimentally determined mean and turbulence parameters are rarely avail- 
able, a very simple core model applicable to simple step-type (slug) profiles was 
developed. Agreement between predicted and experimental mean profiles was gen- 
erally almost as good for calculations made by using this model throughout the core 
region and the transition model for all subsequent regions as predictions made by 
starting from experimental profiles in the transition region. 

The generalized eddy -viscosity model, which was developed in part through 
correlation of turbulence parameters, successfully predicted turbulent shear stress, 
turbulent intensity, and mean velocity profiles for a 0.040-inch-diameter microjet. 
Therefore, successful scaling by the model was demonstrated since data used in its 
development was for jet areas up to 90 000 times as large as the microjet and veloc- 
ities only 1/2 0th as high. 


INTRODUCTION 

There has been considerable interest for many years in turbulence, since most 
flow fields of practical importance are turbulent. Unfortunately, there is no general 
analytical technique available for their prediction. In fact, a group of French scien- 
tists (ref. 1) has objected to the use of a set o f partial differential equations, such as 
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the equations of change, for application to turbulent flows; they feel that since a turbulent 
flow field is in "pure chaos," the instantaneous velocity of a particle of fluid cannot be 
sufficiently regular to satisfy the constraints of partial differential equations such as the 
Navier-Stokes equations. Unfortunately, no substitute for the equations of change has 
been proposed. 

Therefore, the current approach in analysis of turbulent flows is to assume that the 
laminar form of the equations of change apply to instantaneous values of the velocity, den- 
sity, enthalpy, and concentration. These equations are then time-averaged, such that 
cross correlations, that is, Reynolds transport terms, occur in addition to the original 
"laminar -type" terms. Evaluation of these terms requires experimental information con- 
cerning the nature of turbulence, since there are insufficient independent equations to 
specify a given turbulent flow uniquely. The fact that experimental data are required for 
specification of the Reynolds transport terms is the reason that analyses of turbulent 
flows, whatever the approach, must be considered semiempirical. 

Two general techniques that have been used extensively to provide the needed turbu- 
lence input are the eddy-viscosity and the turbulence kinetic energy (TKE) approaches. 

(See refs. 2 to 4.) Both were developed primarily for the prediction of momentum trans- 
port (for example, mean velocity fields) because the velocity field was frequently of pri- 
mary interest. However, as more complex applications arise such as supersonic com- 
bustors and chemical lasers, prediction of mass and energy transport becomes more 
important than momentum transport. At present, the assumption generally made is that 
turbulent Prandtl and Schmidt numbers are constant (generally less than unity, the value 
specified being rather arbitrary) so that the identical approach used for predicting momen- 
tum transport can be used for predicting mass and energy transport as well. 

Unfortunately, the assumption of constant turbulent Prandtl and Schmidt numbers has 
been shown in a number of investigations to be only a rough approximation (for example, 
refs. 5 to 7) so that not until adequate mixing models for mass and energy transport also 
are devised can solution of practical problems involving multispecies turbulent flows be 
obtained with confidence. Of course, generality of an analysis is desirable; however, 
from a practical standpoint, it is not necessary (nor very likely) that all turbulent flows 
will be correlated with a single semiempirical mixing model. As long as flows of a given 
class can be predicted over the complete range of practical interest, useful computations 
can be made. 

Several years ago, the eddy -viscosity approach was selected for further development 
at Bell Aerospace Company, because it appeared to have a number of advantages over the 
newer TKE approach (refs. 8 to 11) for application to practical combustor problems: 
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(1) Successful modeling techniques established for momentum transport could be 
applied directly to the modeling of mass and energy transport, and thereby would make 
the assumption of constant turbulent Schmidt, Prandtl, and Lewis numbers unnecessary. 

(2) Compressible multispecie free shear layer flows were of greatest interest; such 
flows were not successfully treated prior to 1971 by use of the TKE approach. The sim- 
plicity and flexibility of the eddy-viscosity approach for application to practical systems 
such as combustors were judged to be significant. 

(3) Profiles computed downstream of the injection station were less sensitive to 
precise initial conditions. In fact, as demonstrated, the eddy-viscosity approach permits 
calculations to be begun from the injector face by assuming simple (and obviously very 
approximate) step functions for mean velocity, mass faction, and stagnation temperature; 
whereas, the TKE approach requires detailed initial shear stress profiles. 

(4) The eddy -viscosity approach does not require an explicit relationship between 
the shear stress and the TKE. At present there is a controversy as to whether the shear 
stress is directly proportional to the TKE or the square root of the TKE and thus its gen- 
eral applicability is somewhat uncertain. 

(5) Empirical relationships are required for the second and higher order correla- 
tions in the TKE equation. The data used to develop these empirical relationships is 
limited to simple flow conditions, that is, low-speed constant-density flows. Quite pos- 
sibly these relations will not apply in general to more complex practical flow systems. 

(6) The TKE approach had not been demonstrated to predict successfully the details 
of mixing and reacting jet flow. 

(7) Both the TKE and eddy-viscosity techniques use phenomenological approaches, 
and therefore, neither can be expected to describe the detailed physics of the flow. The 
success of predictions and the ease of their use are principal factors to consider in judg- 
ing the relative merits of the approaches. 

SYMBOLS 


D jet diameter 

f ratio of transverse to longitudinal root-mean-square values of velocity 

fluctuations (eq. (2)) 

G radial variation function for eddy viscosity (eq. (9)); also gas 

L characteristic length (eq. (7)) 
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eddy viscosity 


p mean density 

r € Reynolds momentum flux (shear stress) 

Subscripts: 

c evaluated at center line 

e evaluated at r = °° 

j evaluated at r = 0, z = 0 

max maximum value at a given z 

o evaluated at r = o 

u evaluated at velocity half-width 

A bar over a symbol denotes a time-averaged quantity. 

GENERALIZED EDDY-VISCOSITY MODEL FOR COFLOWING STREAMS 

Background 

When the decision was made to pursue an eddy-viscosity approach for turbulent 
mixing analyses, a procedure was selected which, if successful, would circumvent many 
of the shortcomings of previous eddy-viscosity models. Perhaps, most important was 
the decision to model transport of mass, momentum, and energy separately and to con- 
sider the wide range of flow conditions for the model development (refs. 5 and 12) pre- 
sented in table 1. Previous investigators had demonstrated that any select set of data 
could be successfully correlated (refs. 13 to 19) by using various types of models and 
assuming the turbulent Prandtl and Schmidt numbers to be constant. Unfortunately, most 
of these models had rather limited application (refs. 5 and 20). 

The data in table 1 include both very low-speed single-component free jets as well 
as supersonic compressible multispecies coflowing streams for the configuration shown 
in figure 1. Because of the wide range of conditions, it was apparent that some of the 
assumptions made by previous investigators for mathematical simplicity could not be 
valid. For example, the assumption that has been most frequently made for this reason 
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is that no variation in the transverse (radial) direction exists in either the eddy viscosity, 
or alternatively the eddy diffusivity of momentum. (See refs. 13 to 19.) Unfortunately, 
this assumption has no valid basis (ref. 12) so that a general model must be a function of 
radial as well as axial position. 

Another important decision was to apply the computational technique for the numeri- 
cal differentiation of experimental mixing data previously developed (refs. 21 to 23) prior 
to initiating model development. By using this procedure, designated the inverse solution 
technique, mean velocity and concentration profiles are differentiated once in the radial 
direction, each of the remaining terms in an integral form of the shear-layer equations 
evaluated, and the appropriate turbulent transport coefficients solved for as a function of 
position in the flow field. Such information was considered to be essential for successful 
model development. Overall mass and momentum balances were evaluated for consis- 
tency; only high quality data which yielded species mass and momentum balances to within 
20 percent of their average value at each axial station were used for the modeling. 

General guidelines were established for the model development which, if success- 
fully followed, almost certainly would lead to more general eddy-viscosity-type mixing 
models than those previously reported: 

(1) Attempt to model the structure of the turbulence, for example, shear stress and 
turbulent intensities, not merely mean quantities. 

(2) Predictions should exhibit self-preservation in region IV (fig. 1) which are inde- 
pendent of jet initial conditions. 

(3) Consider the lack of attainment of equilibrium of the mean and turbulence param 
eters upstream of the self-preservation region by using anemometry data in the model 
development. Failure to follow this guideline which permits consideration of flow 
"history" was considered an important reason for lack of generality of previous eddy- 
viscosity models. 

(4) Use of the strong points of past models, for example, the eddy-viscosity model 
should yield predictions in agreement with the successful Prandtl-type model for incom- 
pressible free jet flows in the similarity region. 

(5) Judge the merit of the model by comparison of predicted and experimental pro- 
files at each axial station for which data are available. The ability to predict proper pro- 
file shapes was considered to be more important than merely obtaining agreement with 
experimental data along the center line. 

By using these guidelines and the results obtained from the inverse solution tech- 
nique, models for the maximum values of the turbulent shear stress and eddy viscosity 
were devised as well as their variation in the direction transverse to the flow (radial 
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direction). Details of this development were presented in references 5 and 12; a sum- 
mary is presented herein. 


Model Development 

For coflowing streams, in which there is a predominant flow direction, the shear- 
layer equations apply; they are the same as the boundary -layer equations except for 
initial and boundary conditions. The classical definitions of the turbulent shear stress 
applicable to incompressible flows 


Tg = puv = -pR uv u'v' 



( 1 ) 


Equation (1) was used together with the assumption that the effects of density fluctuations 
are negligible, as is usually done in the analysis of free turbulent shear flows. This 
assumption has been shown to be valid for compressible flows (up to Mach 3) in which 
nolecular weight is constant. (See ref. 24.) Empirical relations were developed for vari- 
es parameters in equation (1), that is, for u', v’, and Ry V , which permitted prediction 
f turbulence quantities such as r £ and u'2, as well as mean profiles. 

The data of references 25 to 30 suggested that the ratio of axial to radial turbulent 
ntensity, as a first approximation, may be expressed only as a function of axial position, 
ind that the ratio of radial to tangential turbulent intensity is essentially unity. 


r = f!/2 and * 1 
u’ w' 


( 2 ) 


where 
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The empirical relationship between turbulent kinetic energy and the shear stress 
suggested in reference 11 was used to relate the eddy viscosity to u'^ 
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Therefore, by using equations (2) and (3) and assuming that the region of maximum shear 
occurs at the half-width of the jet r u 


T e _ pu' 1 2 * 8U/3r 
( T e) u (pu’2) (»U/8r) u 


and 


e_ pu' 2 

E “’(pu' 2 ) u 


(5) 


Applying equations (1), (2), and (3) at the jet half- width r u resulted in an expression for 
the maximum value of the correlation coefficient 



-0.15 


1 + 2f 

1/2 


( 6 ) 


The empirical relation for the maximum value of (u' 2 ) 
guidelines and correlated hot-wire anemometry data was (refs. 


that followed the general 
5 and 12) 



= 0.12 



PU - P e U e r dr 


P L 
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(7) 


where 


L = r u t (D - r u )e-°- 115j 

Substituting equations (2), (6), and (7) into equation (1) yields an expression for the maxi- 
mum value of the eddy viscosity 1 

0.018(1 + 2f)(l + e 4 ’ 6U i/ Ue ) J |pu - p e U e |r dr 

% r 2 (8) 


1 The original model developed in reference 12 included an empirical function of 

density ratio. Subsequent analysis suggested that, in general, this function was not nec- 

essary and has been excluded from the model. Interestingly, Brown and Roshko (ref. 34) 

have shown that the mixing rate of planar shear layers also exhibits an independence of 
density ratio. 
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Variation of the eddy viscosity in the transverse (radial) direction was evaluated by 
using equation (5). The empirical expression G^r/r u j used to represent this radial 
variation was 

i_ = P*' 2 = 1.0&-0.15e" 4 - 6r / ru — G / r/r \ (9) 

€u (pu' 2 ) u 1.0 + 0.05(r/r u ) 7 

Equations (8) and (9) define the generalized eddy -viscosity model developed for the 
transition region. This model, together with the assumption that N ScT and Np r 
were constant, was used for predicting the mixing of the free turbulent shear flows 
selected by the conference data selection committee. 

Discussion of Model 

The eddy-viscosity model successfully followed all the general guidelines at least 
partially, and therefore it is unique. Several of its features deserve comment: 

(1) Radial variation of e and pu’ 2 are included in the empirical function 
e / e u = G(r/ r u) defined in equation (9). This function is plotted in figure 2 along with 
constant and variable density data of references 25 to 28 and 31 to 33 which were used 
for its determination. There are several reasons for the scatter of the data: (a) turbu- 

lence quantities have not become self-preserving at the axial stations for which they were 
available and (b) data obtained by using the inverse solution technique scatter even more 
than the anemometry data since numerical differentiation of experimental data is utilized 
in this technique. Consistency of the injected mass and momentum integral balances to 
better than ±10 percent (rather than ±20 percent) is required to reduce the scatter. 

Nevertheless, the composite results in figure 2 strongly indicate that the eddy vis- 
cosity reaches a maximum at some distance from the center line and that it exhibits an 
intermittent-type behavior at the outer edge of the mixing region analogous to the behavior 
of boundary layers (ref. 2). The data also indicate that e/e u is monotonically decreasing 
for r/r u > 1 and that it reaches a maximum at about r/r u ~ 0.6, that is, in the vicinity 
of half-width of the jet. The exact location of the maximum is difficult to pinpoint because 
of the data scatter; however, the maximum of the function defined by equation (9) occurs at 
r/r u = 0.66. No consistent trend was apparent which correlated with density variation 
across the mixing region. Therefore, as a first approximation, the solid curve repre- 
sented by equation (9) was used to represent the average variation of eddy viscosity in the 
transverse direction. This approximation is obviously superior to the assumption that 
e/e u is constant which was frequently made by previous investigators; however, it is 
apparent from figure 2 that a correlation which includes other parameters than r/r n is 
needed to reduce the scatter. 
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(2) Both the eddy viscosity and turbulent intensity u'/Uj are proportional to the 
mass defect. The mass defect has been successfully used by a number of investigators 
to model the turbulent shear stress over a wide range of flow conditions and test geom- 
etries. It was shown by Clauser (ref. 35) to correlate the eddy viscosity in the outer 
region of a two-dimensional turbulent incompressible boundary layer in the absence of 
pressure gradients. Schetz (refs. 18 and 19) extended the mass defect concept to model 
axisymmetric coflowing streams in the transition region. Zelazny (ref. 36) further gen- 
eralized this concept to apply to both jets in a quiescent atmosphere and coflowing 
streams. It should also be noted that the generalized eddy -viscosity model predicts the 
same center-line velocity decay rate in the similarity region as the Prandtl model. 

(3) Recognition was given the fact that turbulent shear stress is not determined 
exclusively by local mean flow properties at a particular location; that is, until the flow 
is in the self-preservation region (fig. 1), prior development or "history" of the flow must 
be considered, for example, initial conditions and wall effects. An attempt was made to 
include this effect in the model empirically by defining a characteristic length L (eq. (7)) 
and a function f (eq. (2)), both of which varied significantly with axial location in the near 
region up to about 100 diameters downstream (z = 100). 

The fact that L = r u (a function of axial position) was a more appropriate charac- 
teristic length than L = D/2 (a constant) was discussed in reference 36. The more 
sophisticated relation of equation (7) was demonstrated to yield superior predictions for 
the data of table 1. Figure 3 is a plot of L for several of the test cases, and shows this 
parameter reaches a minimum value for z ranging from 6 to 10. Since this parameter 
is empirical, and somewhat arbitrarily defined, its further modification may well lead to 
improved results. 

The function f defined in equation (2) varies linearly from 0.5 at z = 0 to 1.0 at 
z = 100. It accounts empirically for the fact that u' and v' are not in equilibrium in 
the near region of a free shear flow. The additional refinement of considering f a func- 
tion of radial position was deemed unnecessary in light of the spread of the data for u'^ 
shown in figure 2. 

The relations L and f both attempt to allow for the fact that the free turbulent 
shear flows do not attain self-preserving profiles for much of the region of practical 
interest. (See fig. 1.) Obviously, these relations are only rough approximations for this 
effect since free turbulent flows do not exhibit a "universal" law by which turbulence 
parameters and mean properties attain self-preserving profiles. Preturbulence levels 
such as those caused by screens and initial profile shapes caused by splitter plates may 
play a significant role in determining the axial location at which this condition is achieved. 

1 2 

2 For constant-density flows, ^- = iil_ (eq. (6)). 

e u „ »2 
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(4) An empirical function, 1 + exp(-4.6Uj/U e ) which yields better correlation 
when U e » Uj was introduced to broaden the range of applicability. However when 
Uj > U e , this function has a negligible effect on the eddy-viscosity model, for example, 
for the axisymmetric test cases the relation ranged from 1.00 to 1.04. On the other hand, 
a significant improvement in predictions was obtained for jets where Uj « U e such as 
those reported by Zawacki and Weinstein (ref. 25). The amplification provided by this 
function, which reaches a maximum value of 2, may account for the increase in the tur- 
bulent momentum transfer caused by the recirculation regions that occur in these "wake- 
like” flows. 


Core Model 

The generalized eddy-viscosity model was developed by using published data for 
coflowing streams. Unfortunately, there were few data available for the core region 
because of the difficulty in obtaining valid measurements in this region prior to the recent 
development of the laser -Doppler technique. Of course, a model applicable to the core is 
of great practical importance since in practical combustors and chemical lasers, for 
example, the ignition and most of the combustion occur within this region. Unfortunately, 
generally neither initial mean profiles nor shear stress profiles are known. The simplest 
type of initial condition for starting calculations would be a "step-type" or "slug" profile 
in which bulk mean quantities were simply used to characterize each stream. Naturally, 
a core mixing model appropriate to step-type initial profiles might not apply when used 
with experimental profiles. However, such a deficiency would not be important for model- 
ing the transition and similarity regions as long as realistic profiles were generated prior 
to reaching the end of the core. An effort to develop a core model is one task of the cur- 
rent Bell Aerospace AFOSR contract. 

As demonstrated in this paper the generalized eddy -viscosity model was quite suc- 
cessful in the transition region for which it was developed. The simplest possible 
approach was to assume a core model with the same functional form; therefore, as a first 
approximation, the transition model was multiplied by a constant factor (less than unity) 
to correct for the overmixing predicted using it in the core. The constant 0.4 proved to 
be satisfactory, and the "core model" which resulted, that is, e core = 0.4et ranS iti On , was 
far more successful than anticipated. Of course, when applying the transition model to 
the core region, the limits of the integral for the mass defect (eq. (8)) were changed from 
0 to 00 (the extent of the mixing zone in the transition region) to rj and ^ (the extent 
of the mixing zone in the core). Success of this model was demonstrated by the very 
reasonable predictions attained when using either experimental or slug profiles at the 
injection station. It suggests that the general functional form of the transition model may 
be appropriate for the core as well; however, further work on core modeling is required 
before any conclusions can be drawn. 
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Additional efforts to model the core region are in progress at Bell Aerospace with 
recent results reported in reference 37. Also reported in reference 37 is a review of 
methods used in the modeling of turbulent axisymmetric coflowing streams and quiescent 
jets. Concepts such as a universal center-line mass fraction decay exponent, the 
Reichardt hypothesis, and a virtual origin were demonstrated generally to be inadequate 
for characterization of these flows. 

COMPARISON OF PREDICTIONS WITH EXPERIMENTAL DATA 

The generalized mixing model was used to predict the flow fields for the twelve axi- 
symmetric jet cases and the two axisymmetric wake cases listed in table 2. The govern- 
ing shear layer equations were solved numerically in the Von Mises coordinates by using 
an explicit finite-difference method. Three types of predictions were made: (1) Step-type 

(slug) initial profiles were assumed at the injection station (x = 0) and the core model used 
until the velocity on the center line of the jet was less than its initial (core) value; the 
transition model was used thereafter. (2) The first experimental profiles reported down- 
stream of the core were used for the initial profiles (x/D designated on curve) and the 
transition model used exclusively. (3) The first experimental core profile was used for 
those cases in which the boundary -layer effects on the splitter plate were not pronounced 
(x/D designated on curve). (Initial velocity profiles that exhibited pronounced boundary- 
layer effects required a programing change to define the extent of the mixing region; this 
change was not deemed to be warranted in the light of the success obtained by using the 
slug profile.) Generally, three sets of curves are presented in each figure and corre- 
spond to each of these types of predictions. Figures are plotted in the manner specified 
by the data selection committee. 

Test Case 6 (Maestrello and McDaid) 

Predictions of the mean velocity for the high-speed quiescent jet are compared with 
the experimental data in figure 4 for test case 6 (ref. 38), results showed good agreement. 
Since no transition region profiles were available, only predictions from the injection sta- 
tion were possible. The predicted center-line velocities were at most 14 percent greater 
than the experimental velocities. This good agreement indicates that the initial slug 
velocity profile was adequate. However, as expected even better agreement was obtained 
when the actual experimental core profile (x/D = 1.0) was used. This result tends to val- 
idate the simple core model for use with realistic initial profiles. 

Test Case 7 (Eggers) 

Velocity data for the supersonic (Mach 2.22) quiescent jet are compared with the 
predictions in figure 5(a) for test case 7 (ref. 32). The agreement achieved appears to be 
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adequate for most engineering applications; however, the model overestimated the mixing 
rate (maximum error was 24 percent). Note, however, that it had underestimated the mix- 
ing rate in the previous case with a maximum error of 14 percent. Again little difference 
in predictions was observed when experimental core profiles were used and thus the 
results of the previous case, were confirmed. Comparison of predicted axial velocity pro- 
files with the data are shown in figure 5(b) for both the slug and experimental initial pro- 
files. Good agreement was obtained at x/r 0 = 8.0 and adequate agreement at x/r Q = 27 
and 99 for each type of prediction. 


Test Case 8 (Heck) 

Predictions of the center -line velocity and total temperature for the high- 
temperature quiescent jet are shown in figures 6(a) and 6(b) for test case 8 (ref. 39). 

Since both momentum and thermal energy transport are significant in this case, it was 
necessary to specify a turbulent Prandtl number. A constant Np r x = 0-70 was 
assumed, since this value appears to be representative of values reported in the literature. 
Calculations were started from the point of injection using the slug and experimental ini- 
tial pro fi les. The former overestimated the mixing rate for both velocity and total tem- 
perature with a maximum error of 20 percent; the latter resulted in somewhat better 
agreement. In addition, when calculations were started from the transition region pre- 
dictions of both velocity and temperature were in excellent agreement with data (maximum 
error was only 8 percent). 


Test Case 9 (Forstall) 

Predictions of the center-line velocity decay for the coflowing air (with 10 percent 
He tracer) mixing with air streams are shown in figure 7(a) for test case 9 (ref. 40). The 
model underestimates the mixing rate with differences between prediction and experiment 
less than 14 percent. Examination of the initial experimental profile and the assumed ini- 
tial slug profile shown in figure 7(b) suggests one reason for the disagreement. The 
assumed momentum flux at the injection station is considerably larger than the actual 
momentum flux because of significant momentum loss to the splitter plate. Therefore, 
mixing was predicted to be slower than actually observed. As a test of the importance of 
this effect, predictions also were made with a slug profile adjusted so that the momentum 
flux obtained from both the experimental and slug profiles would be equal (fig. 7(b)). Pre- 
dicted center-line velocity decay for this case which is also presented in figure 7(a) is in 
excellent agreement with experimental values. In addition, calculations made by starting 
from the transition region show that good agreement was attained between the data and the 
predictions (within 10 percent). 
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Test Case 10 (Chriss) 

Results for the high-speed, subsonic, coflowing streams of hydrogen mixing with 
air are shown in figure 8 for test case 10 (ref. 31). The center-line predictions obtained 
by using both an initial slug profile and experimental core profiles underestimated the 
mixing (29 percent difference for velocity and 57 percent difference for mass fraction). 
No significant difference was observed between these cases. Calculations started from 
the transition region show that slightly better agreement was obtained for the hydrogen 
mass fraction (within 55 percent) but poorer agreement was obtained for the center-line 
velocity (within 45 percent). The first experimental transition velocity profile, obtained 
at x/D = 5.34, is presented in figure 8(c). Good agreement at this station demonstrated 
again the adequacy of the simple core model. However, these results demonstrate that 
the eddy -viscosity model does not include all the complexities required for exact predic- 
tions. The predictions for this case and that of Eggers (test case 12 (ref. 41)) are the 
poorest of the entire set. Nevertheless, results probably still are adequate for many 
engineering purposes. 


Test Case 11 (Eggers and Torrence) 

Predictions for coflowing air (with 1 percent ethylene tracer) and airstreams are 
presented in figure 9 for test case 11 (ref. 7). These results show good agreement between 
predicted and experimental velocities and are obtained by using the initial slug profile 
(within 8 percent) even though data taken at the injector face exhibited boundary -layer 
effects at the splitter plate. Predictions also were made by starting from the transition 
region for this case, and they showed better agreement (within 5.5 percent). 

Test Case 12 (Eggers) 

Results for the Mach 0.89 inner hydrogen jet mixing with a Mach 1.32 outer air- 
stream are presented in figure 10 for test case 12 (ref. 41). A turbulent Prandtl number 
of 0.9 was used in the calculations as suggested in reference 36. Predictions for velocity 
obtained starting both from the core and transition regions (fig. 10(a)) are somewhat low; 
predictions initially were 40 percent too low, but agreement was considerably better at 
downstream stations. The center-line hydrogen mass fraction decay (fig. 10(b)) was 
greatly overestimated (nearly 75 percent too low at x/D = 6). This overmixing is sur- 
prising, since for the similar conditions of Chriss (test case 10) the model predicted 
undermixing, although it gave reasonable agreement overall. Examination of the schlieren 
photographs of these tests suggest that pressure gradients exist in the near region (ref. 41). 
Of course, since the shear layer equations are used in the analysis, transverse pressure 
gradients could not be considered. The omission of pressure-gradient effects may, at 
least in part, account for the poor agreement in figure 10. 
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Test Case 15 (Chevray) 

Predictions for an axisymmetric wake are shown in figure 11(a) for test case 15 
(ref. 42). The model underestimated the mixing rate for this wake. Note, however, that 
the scale w"3/2 is somewhat misleading since it amplifies small discrepancies. For 
example, at x/D = 18.0, the difference between predicted and experimental velocity is 
15 percent; whereas, with respect to the w~3/2 scale, the difference is 70 percent. 
Chevray also has measured the turbulent shear stress and a comparison between data and 
predictions is shown in figure 11(b) at x/D = 12.0. The predicted and experimental shear 
are in reasonably good agreement at this station; however, comparisons at other axial sta- 
tions were not consistently this good. The ability of the model to predict turbulence quan- 
tities is discussed in the next section. 


Test Case 17 (Demetriades) 


Good agreement between experimental and predicted velocities was obtained for the 
compressible wake data presented in figure 12 for test case 17 (ref. 43). The maximum 
disagreement was only 8.0 percent. This result suggests that the model developed by use 
of jet data exclusively is reasonably valid for wake data as well. 


The predictions for the wake flows were made by omitting the term 1.0 + 
exp(-4.6Uj/U e ), from the model since actual wake data were not used in the development 
of this expression (for which Uj = o). Calculations made by including this term improved 
agreement with the Chevray data but agreement between predictions and the Demetriades 
data became poorer. Additional wake data must be evaluated before an empirical velocity 
ratio expression can be validated. 


Optional Test Cases 

Predictions also were made for the five optional test cases (figs. 13 to 17) which 
utilized axisymmetric geometry. Since the previous cases showed that valid results were 
obtained by starting with the initial slug profile, and since such profiles are the simplest 
to use, predictions were made only with these profiles. 

The agreement obtained between the predictions and the experimental data for test 
case 18 (Wygnanski and Fiedler, refs. 29 and 30) was outstanding; that is, similar 
(region in, fig. 1) velocity profiles were predicted. Agreement for test case 19 (Heck, 
ref. 39) was also quite good (within 20 percent) as it had been for the earlier test case 8 
(Heck, ref. 38). The velocity agreement is not very good (within 33 percent) for test 
case 20 (Chriss, ref. 44); however, the hydrogen concentration for this case was somewhat 
better (about 29 percent). As before, agreement was not good for either velocity or con- 
centration for test case 21 (Chriss, ref. 31) (nearly 16 percent for velocity and 100 percent 
for concentration). Agreement for test case 22 (Eggers, ref. 41) also was not good; the 
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model drastically overestimated the mixing rate as it had for test case 12 (Eggers, 
ref. 41). Lack of agreement is not surprising since the data exhibits an anomalous 
behavior in which the center-line jet velocity becomes less than the free-stream velocity, 
evidence that severe pressure gradients may have occurred. For this reason, an analysis 
using shear-layer equations may not be appropriate. 

PREDICTION OF TURBULENCE PARAMETERS 

In the development of the generalized eddy-viscosity model, an effort was made to 
correlate both turbulence and mean quantities. Recently purchased hot wire and laser- 
Doppler anemometers are being used to characterize the flow from small jets in an effort 
to obtain empirical correlations for jet noise. Such data present an excellent opportunity 
to compare turbulence parameters predicted by using the model directly with experimen- 
tal parameters. If successful in these predictions, the model might be used with some 
confidence in the development of an empirical correlation for noise. 

Results of hot-wire measurements made for a 1.02-mm-diameter (0.040-in.) free 
jet by Baker, Moon et al. (ref. 45) are presented in figures 18 to 21; initial jet velocity 
was 213.4 m/sec (700 ft/sec). The data are very consistent since variation in momentum 

balances was less than 2.5 percent from the mean value, that is, \ pU^r dr ranged 

J 0 

from 0.0512 to 0.0533 N (0.01152 to 0.01198 lbf). Predictions were made by using the 
initial slug profile and the core model in the near region followed by the transition model 
as previously described. Agreement is seen to be very adequate for turbulent shear 
stress t £ , axial turbulent intensity u'/Uj, and mean velocity Uj. In the case of r € 
and u'/Uj, agreement was poorest at the initial station (x = 30.48 mm (1.2 in.)) where the 
x-wire is long (1.27 mm (0.050 in.)) relative to the jet diameter. For this reason, char- 
acterization of the near region, that is, core and transition, are in progress using the 
laser-Doppler technique. The agreement between experimental and predicted velocities 
using the model (fig. 20) is outstanding, even at the initial station. The accuracy of pre- 
dictions of mean values (within 2 percent) is significantly better than the shear stress 
(within 20 percent). This result demonstrates that predicted mean values exhibit a degree 
of insensitivity to the inaccuracies introduced in predicting the shear stress. That is, the 
eddy viscosity is used directly for computation of shear stress, but it merely influences a 
coefficient used in the numerical integration procedure to obtain mean velocity. The 
center-line velocity is plotted as a function of axial position in figure 21. The decay 
exponent is correctly given by -1 for both experimental and predicted results as is appro- 
priate for flow in the similarity region. 
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Additional experiments were made for a sonic free jet; although exact static pres- 
sure matching was not attained at the nozzle exit agreement between predicted and exper- 
imental parameters, it was essentially as good as that shown in figures 18 to 20. 

These results enhance our confidence in the applicability of the model for engineer- 
ing calculations. They also showed the model was useful for scaling since data used for 
its development (table 1) were for jet areas up to 90 000 times as large as the 1.02-mm 
(0.040-in.) microjet and velocities l/20th as high. No modification of the model was made 
for the microjet predictions. 

CONCLUSIONS AND RECOMMENDATIONS 

Generalized eddy-viscosity models developed for the similarity, transition, and core 
regions were presented which did creditable jobs of predicting mean velocity profiles (and 
reasonably well for concentration and temperature profiles when appropriate) for most of 
the 14 axisymmetric free turbulent shear flows selected by the data selection committee. 
The core model was demonstrated to be reasonably valid even when slug (step-type) pro- 
files were assumed at the injector. Results established the validity of the eddy -viscosity 
approach for engineering predictions including practical hardware design and its optimiza- 
tion. (For example, see ref. 46.) 

The model also was shown to predict turbulence shear stress and axial turbulence 
intensity as well as mean velocity for jets varying in area by a factor of 90 000 and veloc- 
ities varying twentyfold. Its applicability over a wide range of jet geometries as well as 
flow conditions was thereby established. 

The model has several important features; it includes (1) a transverse (radial) vari- 
ation of the eddy viscosity, (2) the mass defect, (3) allowance for the fact that turbulent 
shear is not dependent exclusively on local mean flow properties (by defining an empirical 
characteristic length, and a function allowing for axial variation in the ratio of the axial to 
radial turbulence intensities), and (4) allowance for the variation in the ratio of jet velocity 
to external stream velocity. 

The assumption of constant turbulent Schmidt and Prandtl numbers was made in 
order to predict flow fields in which mass and energy transport occurred in addition to 
momentum transport. In order to obtain satisfactory predictions, these parameters were 
varied from 0.6 to 0.9. These results clearly demonstrate the need for separate models 
for mass and energy if realistic predictions are to be made without the benefit of prior 
experimental data. 

The critical guideline used in the development of the generalized eddy -viscosity 
model was application of the inverse solution technique to a wide range of valid experimen- 
tal data so that quantitative results were available for determination of the various func- 
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tional relations and constants. The same technique should be applied for the development 
of turbulent mixing models for mass and energy transport, so that the unacceptable 
assumption that the turbulent Schmidt and Prandtl numbers are constant can be relaxed. 

Examination of those workshop cases in which poorest predictions were obtained 
indicates that the ratio (pU)jy/(pU e ) e = 0.6. This result suggests that further improve- 
ment of the eddy-viscosity model may be possible by including an appropriate relation 
containing this ratio. 

Examination of available free turbulent shear flow data suggests that more detailed 
experimental investigations be conducted. These investigations should include the follow- 
ing effects: (1) pressure gradients and pressure levels other than atmospheric, (2) initial 
conditions at the injection station, (3) heavy gas jets exhausting into light gases, for exam- 
ple, oxygen into hydrogen, and (4) jet to free-stream velocity ratios near unity. 
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TABLE 1.- COAXIAL JET MIXING DATA USED IN MODEL DEVELOPMENT 


[External gas stream was air in all cases] 


Investigator 

Case 

number 

Jet gas 

D, mm 

D, in. 

Pj, kg/m3 

Pj, lbra/ft 3 

Uj, m/sec 

Uj, ft/sec 

( z A>j)min 

( Z / D j)max 

(U)j 

(PU)| 

<Pu2). 


M e 

<U) e 

(pu) e 

<Pu 2 >e 

Chriss (ref. 31) 

*1 

»2 

12.7 

0.5 

87.6 

5.47 x 10“ 3 

1005.8 

3300 

2.95 

14.55 

6.300 

0.560 

3.520 

0.79 

0.42 


*2 

H2 

12.7 

.5 

92.7 

5.79 

975.4 

3200 

.49 

14.50 

4.400 

.390 

1.710 

.79 

.60 


*3 

h 2 

12.7 

.5 

87.8 

5.48 

929.6 

3050 

.42 

20.80 

3.800 

.320 

1.210 

.73 

.66 


*4 

h 2 

12.7 

.5 

85.7 

5.35 

731.5 

2400 

.50 

20.80 

3.000 

.240 

.720 

.57 

.67 


*5 

h 2 

12.7 

.5 

82.7 

5.16 

579.1 

1900 

.41 

12.70 

2.400 

.190 

.460 

.44 

.65 


6 

h 2 

12.7 

.5 

80.7 

5.04 

944.9 

3100 

.50 

19.20 

4.600 

.620 

2.850 

.71 

.42 


7 

h 2 

12.7 

.5 

79.3 

4.95 

746.8 

2450 

.50 

16.30 

3.200 

.410 

1.310 

.56 

.49 


8 

«2 

12.7 

.5 

74.0 

4.62 

594.4 

1950 

.51 

12.40 

2.500 

.300 

.750 

.43 

.50 


9 

Air 

12.7 

.5 

1322.4 

82.56 

286.5 

940 

.50 

14.30 

2.400 

3.600 

8.650 

.87 

.71 

Zakkay (ref. 13) 

*1 

H 2 

7.6 

.3 

48.2 

3.01 

603.5 

1980 

10.0 

30.00 

1.460 

.047 

.069 

.51 

1.6 


2 

h 2 

7.6 

.3 

64.1 

4.00 

701.0 

2300 

13.3 

30.00 

1.690 

.072 

.122 

.60 

1.6 


3 

h 2 

7.6 

.3 

76.9 

4.80 

1002.8 

3290 

13.3 

30.00 

2.420 

.124 

.300 

.89 

1.6 


*4 

He 

7.6 

.3 

141.0 

8.80 

454.2 

1490 

13.3 

30.00 

1.100 

.103 

.113 

.51 

1.6 


5 

He 

7.6 

.3 

163.4 

10.2 

691.9 

2270 

13.3 

26.70 

1.670 

.185 

.309 

.82 

1.8 


6 

Ar 

7.6 

.3 

1633.8 

102.0 

219.5 

720 

13.3 

30.00 

.530 

.590 

.312 

.82 

1.6 


7 

At 

7.6 

■3 

2066.3 

129.0 

231.6 

760 

13.3 

26.70 

,560 

.790 

.442 

.89 

1.6 


8 

At 

7.6 

.3 

2354.6 

147.0 

256.0 

840 

16.7 

30.00 

.620 

.976 

.605 

1.00 

1.6 

Eggers (ref. 32) 

*1 

Air 

25.4 

1.0 

2402.7 

150.0 

539.2 

1769 

.0 

75.00 




2.22 

.0 

Eggers and Torrence (ref. 7) 

*1 

Air/ethylene 

24.4 

.96 

1457.6 

91.0 

289.6 

950 

.0 

49.00 

.740 

.640 

.447 

.90 

1.3 

Alpinieri (ref. 47) 

1 

co 2 

50.8 

2.0 

2162.4 

135.0 

95.1 

312 

5.25 

12.50 

.470 

.660 

.310 




*2 

C0 2 

50.8 

2.0 

2226.5 

139.0 

128.6 

422 

5.25 

12.50 

.650 

.950 

.617 




3 

co 2 

50.8 

2.0 

2114.4 

132.0 

154.2 

506 

5.25 

12.50 

.780 

1.170 

.912 



Zawacki and Weinstein (ref. 25) 

*1 

Air 

18.14 

.714 

1217.4 

76.0 

.366 

1.20 

.0 

11.30 

.025 

.025 

.001 




*2 

Air 

18.14 

.714 

1217.4 

76.0 

.503 

1.65 

.0 

11.30 

.035 

.036 

.001 




3 

Air 

18.14 

.714 

1217.4 

76.0 

.920 

3.02 

.0 

14.00 

.063 

.063 

.004 




*4 

Air 

18.14 

.714 

1217.4 

76.0 

1.829 

6.00 

.0 

14.00 

.125 

.125 

.063 

i 



*5 

Air 

18.14 

.714 

1217.4 

76.0 

4.328 

14.2 

.0 

21.00 

.294 

.294 

.089 




6 

Air 

18.14 

.714 

1217.4 

76.0 

14.630 

48.0 

1.40 

14.00 

1.000 

1.000 

1.000 




7 

Freon 12 

18.14 

.714 

4869.5 

304.0 

1.207 

3.96 

.35 

14.00 

.086 

.344 

.030 




*8 

Freon 12 

18.14 

.714 

4837.4 

302.0 

2.621 

8.60 

.35 

21.60 

.179 

.714 

.127 




*9 

Freon 12 

18.14 

.714 

4869.5 

304.0 

1.341 

4.4 

.35 

14.00 

.185 

.740 

.137 




*10 

Freon 12 

18.14 

.714 

4837.4 

302.0 

.564 

1.85 

.35 

14.00 

.132 

.526 

.070 



Sami (refs. 26 to 28) 

*1 

Air 

304.8 

12.0 

1217.4 

76.0 

10.668 

35.0 

1.00 

10.00 






Wygnanski (refs. 29 and 30) 

tl 

Air 

25.5 

1.004 

1217.4 

76.0 

J 

57.912 

190.0 

5.00 

100.00 







*(1) Momentum integral balances were within 4 percent of their average value at any axial station. 

(2) Injected mass integral balances were within 20 percent of their average value at any axial station and were within 13 percent of their average value at 
80 percent of the axial stations. 

t Detailed profiles not available. 
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Figure 1.- Schematic of coaxial turbulent jet and definition of mixing regions. 
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*See Table 1 



Key 

Sym Investigator/Case*/ i 

O Sami/1/6 
□ Eggers/1/12.6 

O Chriss/1/10.3 
O Chriss/2/12.5 

V Zawacki and Weinstein/ 9/5.6 
A Zawacki and Weinstein/9/5.6 
A Zawacki and Weinstein/10/8.4 

• Sami/1/10 
«T Sami/1/10 
^ Gibson/1/50 

m Zawacki and Weinstein/1/1 1.2 

♦ Zawacki and Weinstein/2/1 1.2 

# Zawacki and Weinstein/4/11.2 

V Zawacki and Weinstein/5/1 1.2- 

Solid Symbols obtained from 
Anemometry Data for Constant 
Density Flows. 

Open Symbols obtained by 
Differentiation of Mean Velocity 
Data for both Constant and Variable 
Density Flows 


Figure 2.- Radial variation of e ^e u . 













° r c = 5.46 cm (2.15 in.) Test Case 8, Heck (Ref 39) 

^ r o “ •318 cm (0.125 in.) Test Case 9, Forstall (Ref 40) 

© r c = 3.10 cm (1.22 in.) Test Case 6, Maestrello and 

McDaid (Ref 38) 

a r 0 = 1.28 cm (0.503 in.) Test Case 7, Eggers (Ref 32) 

A r 0 = 0.635 cm (0.250 in.) Test Case 10, Chriss (Ref 31) 


10 20 30 40 50 

x/D 


Figure 3.- Predicted characteristic length (eq. (7)). 
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u 0 = 211.2 m/sec (693 ft/sec) 
r Q = 3.10 cm (1.22 in.) 


Figure 4.- Predicted and experimental center-line velocity for test case 6 
(Maestrello and McDaid (ref. 38)). 
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(b) Axial velocity against distance from center line. 
Figure 5.- Concluded. 


304 



(b) Center-line total temperature. 

Figure 6.- Comparison of predicted and experimental velocities and total temperatures 

for test case 8 (Heck (ref. 39)). 






O Exp., Forstall (Ref 40) 

Predicted from x/D = 0 using "Slug" Profile 

Predicted from x/D = 8.1 using Exp Profile 

— Predicted from "Slug" Profile with Momentum 
Matched to Exp. Profile (See Figure 7b) 


u Q = 36.57 mps (120 fps) 
D = 0.635 cm (0.25 in.) 


(a) Predicted and experimental center-line velocity 
Figure 7.- Velocity and slug profiles for test case 9 (Forstall (ref. 40)) 
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O Exp., Chriss (Ref 31) 



(a) Predicted and experimental center-line velocity. 

Figure 8.- Velocities and H 2 mass fraction for test case 10 (Chriss (ref. 31)). 
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(b) Predicted and experimental center-line H2 mass fraction 
Figure 8.- Continued. 
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(c) Predicted and experimental axial velocity as a function of distance from center 

line at x/D = 5.34. 


Figure 8.- Concluded. 
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O Exp., Eggers and Torrence (Ref 1) 
Predicted 


u e * 390.0 mps (1279.5 fps ) 
D = 2.44 cm (0.96 in.) 



Figure 9.- Predicted and experimental center-line velocity for test case 11 

(Eggers and Torrence (ref. 7)). 
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U 0 = 1074.1 mps (3524.0 fps) 
D = 1.16 cm (0.456 in.) 


O Exp., Eggers (Ref 41) 
Predicted 



(a) Predicted and experimental center-line velocity. 

Figure 10.- Velocities and H 2 mass fraction for test case 12 (Eggers (ref. 41)). 






O Exp., Eggers (Ref 41) 


D = 1.16 cm (0.456 in.) 

^Sc,T “ 9-9 Predicted 



Figure 10.- Concluded. 
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6 8 10 12 14 16 18 

x/D 


(a) Predicted and experimental center-line velocity. (1 - u/u e )~ 3 / 2 as a function of x/D. 
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(b) Predicted and experimental shear stress profile at x/D = 12.0. 

Figure 11.- Velocities and shear stress profiles for test case 15 (Chevray (ref. 42)). 
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U e * 618.1 mps (2029 fps) 
D = 0.396 cm (0.156 in.) 


O Exp., Demetriades (Ref 43) 
— Predicted 



x/D 

Figure 12.- Predicted and experimental center -line velocity for test case 17 

(Demetriades (ref. 43)). 



0 Exp., Wygnanski and Fiedler (Refs, 29 and 30) 
— Predicted Using Slug Profile 



o 2.0 4.0 6.0 8.0 10.0 12.0 14.0 16.0 18.0 20J 

x/D 

Figure 14.- Predicted and experimental center-line velocity for test case 19 

(Heck (ref. 39)). 
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0 Exp., Chriss and Paulk (Ref 44) 

Predicted Using Slug Profile 

N ScT =0.65 


0 2.0 4.0 6.0 8.0 10.0 12.0 14.0 16.0 18.0 2 

x/D 

(a) Predicted and experimental center-line velocity. 


© Exp., Chriss and Paulk (Del 44) 
— Predicted Using Slug Profile 
T = 0.65 


0 2.0 4.0 6.0 8.0 10.0 12.0 14.0 16.0 18.0 2 

x/D 

(b) Predicted and experimental center-line H 2 mass fraction. 

Figure 15.- Velocities and mass fraction for test case 20 (Chriss and Paulk (ref. 44)) 





(a) Predicted and experimental center-line velocity. 

Figure 16.- Velocities and mass fraction for test case 21 (Chriss (ref. 31)). 
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Figure 17.- Predicted and experimental center-line velocity for test case 22 

(Eggers (ref. 41)). 
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poise ft -sec 2 
sen 40 r 


T « 1.5 


x/D = 30 


u„ = 229 mps (700 fps) 

D = 0.1015 cm (0.04 in.) 

® x/D = 30 ) 

A x/D = 50 > EXP., BAKER, MOON, et al (Ref 45) 
□ x/D = 70 J 


Predicted 


x/D = 50 
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□ \D 


x/D = 70 


OL 0 1 1 lA 1 N. iN. i i 

0 5 10 15 20 25 30 

r/r Q> Distance from Centerline 

Figure 18.- Shear stress profiles at x/D = 30, 50, and 70 for data of Baker, 

Moon, et al. (ref. 45). 



Uj = 229 mps (700 fps) 

D = 0.1015 cm (0.04 in.) 
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O = x/D = 30 
A= x/D = 50 
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| EXP -r 


Baker, Moon, et al (Ref 45) 



Figure 19. - Predicted and experimental axial turbulent intensity profiles at 

x/D = 30, 50, and 70. 
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Figure 20.- Predicted and experimental axial velocity profiles at x/D = 30, 50, and 70. 
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Figure 21.- predicted and experimental center-line velocity. 
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DISCUSSION 


S. Corrsin : This is a trivial remark, but on the figure showing eddy viscosity e versus 
radial distance r you have a finite slope of e at the center line of the axis of the jet. 

I wonder whether it would be better to have a horizontal slope there. 

J. H. Morgenthaler : Of course, the center line is an axis of symmetry; however, because 
(Su/Sr)^ = 0, it is not possible (or necessary) to define (9e/9r)^. 

I. E. Alber : Also, with respect to your eddy viscosity model, I noticed that you had the 
eddy viscosity e proportional to u'2. Nominally, from dimensional analysis, one would 
take e proportional to u\ Is there any reason that you chose that over the more con- 
ventional scheme? 

J. H. Morgenthaler : Steve, would you like to answer that one? 

S. W. Zelazny : The proportionality you are referring to (eq. (5)) was obtained by assuming 
(1) the shear stress is directly proportional to the turbulence kinetic energy, an assump- 
tion used in a number of the Workshop papers using the turbulence kinetic energy 
approach, (2) the ratio of the transverse turbulence intensity to the axial turbulence inten- 
sity is independent of radial position, and (3) the shear correlation coefficient is directly 
proportional to the partial derivative of the axial velocity with respect to the radial coor- 
dinate, that is, 9U/9r. The relation, e proportional to u'2, also may be obtained from 
dimensional analysis if it is assumed that the eddy viscosity is functionally dependent on 
density u'2 and the maximum value of 9U/9r at a given axial station. 

J. H. Morgenthaler : In other words, we did not pull the relation out of our hat. 

M. V. Morkovin : Concerning the same point, how do you use equation (5)? Do you use e 
just to evaluate u'2? It is not used in the development of other relationships, is it? 

S. W. Zelazny : Equation (5) was used to develop the empirical expression, G(r/r u ) of 
equation (9) describing the radial variation of eddy viscosity and turbulence intensity. It 
showed that both eddy viscosity "data" and turbulence intensity data could be used to 
obtain an approximation for e/e u and (pu'2)/(pu'2) u consistent with the assumptions 
listed in my reply to Dr. Alber. It is not true that e is just used to evaluate u'2 since 
e is essential for prediction of the mean values. 

J. H. Morgenthaler : In figure 2 of the paper, there are open symbols and closed symbols 
which didn't show up very well on the slide, but as you can see, some points were obtained 
by direct differentiation of the mean data, and some were obtained directly from hot-wire 
turbulence measurements. 
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M. V. Morkovin : But I’m saying, that all you are trying to do is to fit some transverse 
variation of e. In your computations as you are marching in the axial direction, the 
effect of transverse variation does not influence predictions since the function G(r/r u ) 
is independent of z. That is, the transverse variation only affects prediction of the tur- 
bulence quantities, doesn’t it? 

J. H. Morgenthaler : The transverse variation of e has a great influence on predictions, 
as is illustrated in the figure. The local value of e greatly influences the rate of 
momentum transport as minor perturbations in e have demonstrated during the devel- 
opment of the model. It is true, however, that the computation of mean quantities is not 
influenced by the prediction of turbulence intensity. In other words, the model may pre- 
dict valid mean profiles but be considerably poorer in its predictions of turbulence 
quantities. 

T, Cebeci : Well, let us have one more question before we take our coffee break. 

S. C. Lee : I am looking at your figure 2 right now. I have a question related to this fig- 
ure. It looks from the figure that the function G(r/r u ) might be considered to be con- 
stant between r/r u = 0 and r/r u = 1.5. What effect would this assumption have on the 
predictions? 

J. H. Morgenthaler : We believe the radial variation to be important based on our model- 
ing experience. We have run a case which shows that the length of the transition region 
is significantly influenced as are mean quantities (see the figure). For example, the 
velocity at the t at z = 50 D (2.0 in.) was 92.0 ft/sec for G(r/r u ) defined by equa- 
tion (9), but was 64.3 ft/sec when the assumption was made that G = 1. 
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45.72 r 


150 


Uj = 213.36 m/sec (700 ft/sec) 
D = 0.10 cm (0.04 in) 



Figure 1.- Predicted velocity profiles at z/D = 50 showing effect of 
including radial variation in eddy viscosity model. 
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